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Abstract 

We consider single-channel transmission through a double quantum dot system that consists of 
two single dots coupled by a wire of finite length L. In order to explain the numerically obtained 
results for a realistic double dot system we explore a simple model. It consists, as the realistic 
system, of two dots connected by a wire of length L. However, each of the two single dots is 
characterized by a few energy levels only, and the wire is assumed to have only one level whose 
energy depends on the length L. The transmission is described by using S matrix theory. The 
model explains in particular the splitting of the resonant transmission peaks and the origin of the 
transmission zeros. The latter are independent of the length of the wire. When the transmission 
zeros of the single dots are of first order and both single dots are identical, those of the double 
dot are of second order. First-order transmission zeros cause phase jumps of the transmission 
amplitude by vr, while there are no phase jumps related to second-order transmission zeros. In 
this latter case, a phase jump appears due to a resonance state whose decay width vanishes when 
crossing the energy of the transmission zero. 



I. INTRODUCTION 



Advances in nanotechnology have made it possible to fabricate quantum dots (QDs) and 
to study their transport properties. Single quantum dots can be considered as artificial 
atoms if the energy levels can be resolved. Two or more QDs can be coupled to form 
an artificial molecule, in which the electrons are shared by different sites. Recently, much 
interest is devoted to the study of the properties of these coupled QDs or artificial molecules 
Q]. One of the most important challenges is to understand the basic properties of coupled 
QDs since they display the simplest structures of quantum-computing devices that can be 
controlled by means of external parameters. An example is the interdot coupling that can 
be tuned by external parameters far out of the regimes known for natural molecules. The 
round-state interatomic distance is dictated, in natural molecules, by the nature of bonding 
^ while such a restriction does not exist in artificial molecules. Most of the work has focused 
on the situation where the coupled QDs are in the tunneling regime so that the physics is 
dominated by Coulomb blockade effects 0, 0] . It was observed a pronounced anticrossing 
resonance scenario (called also avoided level crossing scenario) instead of a simple crossing 
one. Recently Rushforth et al. studied the movement of electrons between the dots. 
Application of a voltage onto the gate between the dots allows to vary the length as well as 
the width of the wire connecting the dots. The study revealed evidence that electrons move 
between the dots via excited states of either the single dots or the double dot molecule. 

In the present paper, we consider the ballistic transport through a double QD in the 
regime where Coulomb blockade effects can be neglected and where the transmission is 
resonant. The two single QDs are connected by a wire of finite length. This gives the 
possibility to vary the length or the width of the wire and to study the transmission through 
the double QD as a function of the wire's size and energy. When the length L of the wire 
is much larger than its width then the wire has at least one eigenmode with the energy 
e oc This mode appears additionally to the eigenstates e oc of the two single 

QDs (where R is the radius of the dot) which may be equal to one another or different 
from one another. Such a double QD system allows therefore to investigate, among others, 
the quantum mechanical problem of the coupling of two identical quantum systems that is 
provided by a third quantum mechanical subsystem with an own energy spectrum. This 
system is the analogue of a molecule with hydrogen bonds. 
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Some years ago, the phase of the transmission amphtude has been measured in a double- 
sht interference experiment j^. The results showed phase jumps by vr between resonances 
which raised intensive theoretical work for an explanation 13,0,0,^, U, 3l- Most of these 
calculations associate the sharp phase drops with the occurrence of transmission zeros and 
relate them to the interference zeros of Fano resonances. In ^^l, it was shown however that 
the existence of a transmission zero is, indeed, a necessary condition for the phase jump but 



121, on the destructive 



not a sufficient one. The sharp phase change bases, according to 
interference between neighboring resonance states and thus differs from the mechanism based 
on the Fano interference picture. Destructive interferences between neighbored resonances 
reconsidered also .nflfl. 

The Fano resonance phenomenon characterizes the interference between a single reso- 
nanee with a .e.ative.v s,„ootK Q. TKe interference „ in the reginre 

of overlapping resonances are, however, much more complicated than those in the regime 
of isolated resonances. This has been demonstrated, e.g., in an experimental study of the 
conductance through a quantum dot in an Aharonov-Bohm interferometer uM and in a 

n n 

theoretical study [1^. These results are another hint to the conclusion, drawn in that 



the sharp phase changes observed in p] are the result of processes being different from the 
simple Fano interference picture. 

Here, we study the transmission properties of the double QD system when one lead is 
attached to the first single QD, another one to the second single QD, and both single QDs 
are connected by a wire of finite length L. Also in this system, transmission zeros appear. 
The special situation of a double QD is such that the transmission zeros of the whole system 
are determined by the zeros in the transmission through the single QDs. Since phase jumps 
in the transmission amplitudes are related to the transmission zeros, the mechanism of their 
appearance in a double QD system is expected to be different from that based on the simple 
picture of Fano resonances. 

We will show in the present paper that transmission zeros of first order of a single 2d 
QD cause phase jumps by vr in the amplitude of the transmission through this single QD. It 
depends on the spectral properties of both single QDs and on the manner they are connected 
to the wire and to the leads, whether or not the transmission zeros of the single QDs cause 
phase jumps of the transmission amplitudes of the double QD system. When the single QDs 
remain true 2d-systems in the double QD (as shown in Fig. |^ and have different energy 
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spectra, each transmission zero of each single QD causes a transmission zero of the same 
type in the double QD system. When the spectra of both single QDs are however equal 
(and their connection to the wire and the leads is the same as above), the corresponding 
transmission zeros of the double QD system are of second order. They give rise to two 
phase jumps, each by vr, that compensate each other. When a resonance state crosses this 
transmission zero at a certain length of the wire, its decay width vanishes and a phase jump 
appears now due to the extremely narrow resonance. The transmission zeros do not depend 
on the length of the wire when there is only one channel for the propagation of the mode in 
the wire and in the leads. 

In Sect. II, we present some numerical results for the transmission through a double QD 
that consists of two single QDs connected by a wire. The main features of the transmission 
are represented as a function of the length L of the wire inside the double QD: the trans- 
mission zeros are independent of L while many transmission peaks show some periodicity as 
a function of L. In the following sections, we study the transmission through a double QD 
in detail by using a simple model with only a few levels in both single QDs and one level 
in the wire. The description is based on the S matrix theory for the transmission through 
quantum dots (l^. In Sect. Ill, the formalism is derived and some typical numerical results 
are given and discussed. The spectral properties of the double QD are characterized by 
the (complex) eigenvalues of the effective Hamiltonian that describes the double QD when 
opened by attaching the two leads to it. This effective Hamiltonian is non-Hermitian, and 
its eigenvalues provide both the positions in energy and the decay widths of the states. The 
appearance of transmission zeros is determined by the spectral properties of the two single 
QDs. 

The relation between transmission zeros and phase jumps is discussed in Section IV. 
Most results are obtained for the case that the two single QDs keep their 2d structure when 
connected to the leads and to the wire. The results show very clearly that the mechanism 
differs from that based on the simple Fano interference picture: the spectral properties of 
the states of the double QD system depend strongly on the length of the wire, while the 
transmission zeros do not depend at all on the length. The transmission zeros are of first 
order when the energy spectra of the two single QDs are different from one another, and of 
second order when the spectra of the two single QDs are identical. With the last section, 
we conclude the present study by mentioning some results for more complicated double QD 
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systems. In particular we consider the case that the widths of the leads and of the wire are 
different from one another. When the wire's width exceeds the lead's one, a few channels 
can propagate through the wire and the transmission zeros depend on the wire's length. 



II. SINGLE CHANNEL TRANSMISSION THROUGH DOUBLE DOTS 

Let us consider first a double QD consisting of two identical single QDs, for example, two 
circular dots, that are connected by a wire of the length L. For simplicity we choose the 



width of the wire d equal to the width of the leads, see Fig. 
through such a system is given by Klimenko and Onipko 
probability of the transmission for the single-channel case is 



'a). The theory of transmission 
1. If the dots are identical, the 



(1) 



Ti2 + 4(1 -ri)sin2(0 + A;L) 

where Ti = ti is the complex amplitude of the transmission through the single 

QD, (j) = arg{ti), and k is the wave number related to the energy of the single- channel 
transmission as 

E = k^ + 7t\ (2) 

All values are dimensionless via the characteristic energy ^^,^2 and the width d of the wire 
and the leads. In the expression (^, evanescent modes are ignored whose wave vectors are 
imaginary. That means, (0) can be used for kL ^ 1. The results shown in Fig. ^ are 
obtained from numerical calculations with A; > 3. Here, formula is applicable already for 
L > 1. Similar results for a system of two single QDs with different shape were obtained 
numerically by Pichugin [isj ]. 

The transmission probability demonstrates a few interesting features that are shown in 
Fig- Q (b) - (d). The first is some periodical dependence of the transmission peaks on 
the wire length and on the energy. This dependence can be seen immediately in formula 
(^, and we will not discuss it further. The second feature are the transmission zeros. 
They do not at all depend on the wire length, see Fig. Q (d) where the transmission 
probability T is shown in logarithmic scale. Also this feature follows directly from (Q): the 
zeros of the transmission probability Ti through the single QD lead directly to zeros of the 
total transmission probability T. The third feature concerns the peaks of the transmission 
probability. Fig. ^ demonstrates crossings and anticrossings of resonant peaks as observed 
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FIG. 1: (a) The shape of the double QD that consists of two identical circular dots and a wire 
connecting the two dots. Numerical size of the double QD: radius i? = 40 and width of the wire 
and leads d = 8. (b) The probability Ti{E) for the transmission through one of the single circular 
dots versus energy, (c) The probability T[E, L) for transmission through the double QD shown in 
(a) versus energy and length of the wire, (d) The same as (c) but ln{T{E, L)) in order to show 
clearly the zeros for the transmission through the double QD. 

in experiments . Some of them are independent of the length L, as the transmission 
zeros, while other ones are dependent on L. The latter ones cross the L independent peaks 
as well as the zeros, but the behaviour in the vicinity of the crossing with a zero is different 
from that in the vicinity of the crossing with a maximum. 

In the following, we will consider these features in detail. For this purpose, we use the 
periodicity (first feature) of the transmission picture that allows us to restrict the investiga- 
tion to the transmission properties of a simple model with only a few states. The study is 
based on the S-matrix theory. 
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III. S-MATRIX DESCRIPTION OF THE TRANSMISSION THROUGH DOUBLE 
DOTS 



A. Closed system consisting of two single dots connected by a wire 



The Hamiltonian of the double QD shown in Fig. ^ (a) consists of three parts: two parts 
describe the two single QDs and a third one is related to the wire. The Hamiltonian of the 
two single QDs is formulated in standard way, 



Nl Nr 

Hd = Hl + Hr = £nML){nL\ + enR\nR){nR\ 

nL=l riR=l 



(3) 



where the indices L, R stand, respectively, for the left and right single QD with the energies 
£nL,£nR and the Hilbert dimensions Nl,Nr. The wire is the third independent quantum 
mechanical subsystem described by the Hamiltonian 



H,, 



nu,=l 



{L)\n^){n 



(4) 



We assume that the eigenenergies e„^(L) of the wire depend only on its length L. This offers 
the possibility for the energies of the wire to cross the eigenenergies of the single QDs. We 
assume further that the wire is coupled to, respectively, the left and right single QD via the 
matrices Ul, Ur of rank A'^^: x Nr x dw Then the total Hamiltonian has the following 
matrix form 



Hl Ul 
Ut Ur 



\ 



(5) 



V Ur Hr) 

The Hamiltonian (jS)) differs from those used in the literature for the description 

of a double QD of similar shape by taking explicitly into account the third part (j3]) for the 
wire. 

For the simplest case A^^^ = A^^ = Nr = 1 and equal single QDs, the total Hamiltonian 
takes the following form 

^ €i U ^ 

Hb= u t{L) u . (6) 

U El I 



V 
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The eigenvalues of this Hamiltonian are 

ei + e{L) 



El,. 



7] 



TV, E2 = £i, 



(7) 
(8) 



and the eigenstates read 

/ 

1 

|1) = 



277(77 + Ae) 



-u 



T] + Ae 

—u 



/ 1 \ 




|3) 



2ri{ri - Ae) 



u 

1] — Ae 
u 



(9) 



It is remarkable that one of the eigenenergies of the total system coincides with the energy 
ei of the single QD. This fact remains also in the more general cases with higher dimensions. 

Let us consider the case of two identical single QDs with the number N of states that are 
coupled via the wire. Then it follows Ul = Ur, if the wire is a straight one. The determinant 
which defines the eigenvalues of the total Hamiltonian is of rank 2N + iV^, 



Ei-E Q 
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Un 
U12 






E'y — E • • • U' 



21 



Ul2 
U22 



... 
... 








U21 
U22 






■■■ ei{L)-E ■■■ U12 Un 

■■■ e2{L)-E ■■■ U12 Un 



U21 
Un 



U22 
Uu 



■ ■ ■ E2 
... 



E 
El 



E 



0. 



(10) 



When E = En, n = 1, . . . ,d, two lines in (fTTHl coincide. Therefore, among the 2N + 
Ny^ eigenenergies of the Hamiltonian (jSl), there are N eigenvalues which are equal to the 
eigenenergies of the single dots. Numerical examples for the eigenvalues of Hb with = 1 
and 2 are given by the solid lines in Fig. El 



B. S-matrix for the transmission through double dots: single QDs with one state. 

The knowledge of the eigenstates of the closed quantum system allows us to formulate 
the S-matrix and the effective Hamiltonian in the manner described in Q]. Let us consider 
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w 



u 




u 




FIG. 2: Two single state QDs are connected to the wire w with the couphng constants u and to 
the reservoirs with the couphng constants v. 

first the most simple case of the Hamiltonian ©. For the transmission through this system 
we imply two leads coupled to the reservoirs with the strength v, as shown in Fig. |21 
The coupling matrix has the following general form 16|, |20| 

^ = E E fdEVm{E,C)\E,C){m\ + H.C. (11) 

m C=L,R 

where \m) are the eigenstates of the closed system given in the present case by 0, and C 
enumerates the reservoirs with the states \E,C) normalized by 

{E,C\E',C') = 5{E-E')5c,c'. 

Obviously, 

VUE,C) = {E,C\V\m) = ^ / dxMx)i^Uj) {^\V\j) (12) 

where j runs over the double dot system, x spans over the C-th reservoir, ipc{x) are the 
eigenf unctions of the reservoirs, and ipmU) are the eigenf unctions of the closed double QD 
system. In our model we choose the couplings, that are shown in Fig. |21by solid lines, to 
be fixed at some points: at j = 1, 3 of the left and the right single QD and at some points 
xc which belong to the reservoirs. Moreover, we describe for simplicity the reservoirs as 
semi infinite one-dimensional wires in tight-binding approach 

Q- As in Q, 

we take the 

connection points of the coupling to the reservoirs at the edges of the one- dimensional leads. 
Then the matrix elements take the following form 



Vm{E,L) = VlijE,Li^L)lpm{j = 1) 



VmiE,R) = vi)E,L{xR)i)m{3 = 3) = f 1 



/sin k 



^m(3), 



(13) 



where k is the wave vector related to energy hy E = —2 cos A; + 2. For continual case the 
last equality is simply E ^ k"^. The effective Hamiltonian can be written as 

1 



16 



2Q| 



H, 



cS 



Hb+ J2 ^1 



C=L,R E+ - He 



Vc 



(14) 



where He is the Hamiltonian of the reservoir C and E'^ = E + iO. Substituting (|T!?|l into 
(fT^ we obtain for the matrix elements of the effective Hamiltonian 

^,VUE',C)V^{E',C) 



1 r 

{m\HeG\n) = EmSmn+ ^ 

^ r r, 27r Jn 



dE' 



C=L,R 



E + iO-E' 



= EmSmn ' (^^^ (1)^„ (1) + ^^(3)^„(3))e^^ (15) 

where the states ipn{j) are given in and the indices j = 1,3 mean, respectively, the left 
single QD (j = 1) and the right single QD (j = 3). Substituting Q into (fT3|) we obtain 



H 



eff 



E, 







V 



El- v'^e^^ 









E' 



(16) 



I y-u e 

3 »7(r)— Ae) / 



Next we calculate the (complex) eigenvalues of the effective Hamiltonian (|16|) that are related 
to the poles of the S matrix. The result is 



Z2 



El + e(L) - v^e''' 
^1,3 = -i: T 



\ 



+ 2^2. 



(17) 



For u = 0, the eigenvalues are zi^s = Ei — f^e*'^, Z2 = ^(L). The eigenvalue Z2 = e(L) means 
that the wire has no other connection to the reservoirs as that via the single QDs. 

In order to calculate the S-matrix we need the eigenstates of the effective Hamiltonian 

Q 

i^efflA) = za|A), (18) 

where (A| = |A)'^, A = 1,2,3, c means transposition, and (A|A') = (5a,a' is the biorthogonality 



relation for the eigenfunctions of the non-Hermitian Hamiltonian H^s 
from (fT^ (for u ^ 0): 



2l|. Then if follows 



|1) 



a 




|2) 



1 



|3) 



/ b \ 




(19) 
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with 

a = - , ^ , b=\m^ (20) 

and 

The knowledge of the eigenstates (fTTHl of the effective Hamiltonian allows us to write the 
amplitude for the transmission through the double QD in simple form 



3, 



,^_,„j.(mmm. (22) 

A E - Zx 

The transmission probability is T = Substituting (fTSj) . (fT^ and correspondingly (jH)) 
into the matrix elements (L|y|A) and (A|\^|-R) we obtain 



{L\V\2) = j:{E,L\V\m){m\2) = lf-^, 



i2\V\R) = Y.mm){m\V\E,L) = -lf-^., 



{L\V\1) = il\V\R) = vf-^iMl)a + Mm, 



{L\V\3) = {3\V\R) = v^'-^{Ml)b-Ml)a), (23) 

where the eigenstates ipmij) are given in (jUj). Substituting finally into we obtain 
the transmission amplitude t. 

The typical behavior of the transmission probability T = \t{E, L)p versus energy E and 
length L is shown in Fig. IHl Here the energies Sk of the two single QDs and the energy e(L) 
of the wire are shown by dashed lines while the eigenenergies ((Zj) of the double QD system 
are shown by solid lines. Since the eigenenergy E2 of the double QD system coincides with 
the energy ei of the single QDs, the last one is not shown in Fig. |21(b). The eigenenergy of 
the wire is assumed to depend on the length L according to e(L) = —1/2 — L/5. This linear 
dependence is not decisive for the following discussion, see the Concluding Remarks. 

We underline that the results presented in Fig. IHl follow from a simple model which 
describes the double QD system by two single-state dots connected by a wire. The wire 
is characterized by the only energy e(L). When the coupling of the double QD system 
to the reservoir is relatively weak [meaning that the ratio v/u is small as in Fig. El (a, 
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-ij6 -i.i -12 -1 -as -as -15 -1 -as 

E E' 




riSy --1 -as -15 -1 -as 

B E. 



FIG. 3: The transmission probability T through the double QD shown in Fig. [3 (a) £i = —1; v = 
0.3, u = 0.2 and L = 4. (b) The same as (a) but the length L is not fixed, (c) The same as (b) 
but £1 = -1, v = 0.6, u = 0.2. (d) The same as (b) but e{ = -1.2, ef = -0.8, v = 0.3, u = 0.2. 
The eigenenergies Ek{L), k = 1,2,3, of Hb are shown by full lines while the energies ei and 
e(L) = —1/2 — L/5 are given by dashed lines. 

b, d)], the transmission probability follows the eigenenergies of the closed double QD, and 
we have resonant transmission. We can now compare Fig. 01 (b) with Fig. Q and we see 
that the simple model, basic of Fig. |21 (a), is able to explain the peaks of the transmission 
probability. That means, the main features of the transmission picture in Fig. ^ are related 
to the crossing of the levels of the single QD with the energy level of the wire. As shown in 
Fig. ini(b), (c) and (d) different scenarios of the resonant transmission peaks can be realized 
from the crossing to the anticrossing behavior in dependence on the parameters, basically 
on the coupling strength between the system and the attached leads and on the coupling 
strength between the dots and the wire. 

We can further learn from Fig. El that the transmission through a system consisting of 
two dots that are connected by a wire, is characterized by the ratio v/u. Here, v is the 
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coupling strength of each single QD to the corresponding lead and u is the internal bonding 
of the two single QDs. Thus, v/u is the ratio between external and internal coupling of 
the states of the double QD via, respectively, the reservoir and the wire. In Fig. |21 (c), the 
coupling V between the double QD and the reservoirs exceeds essentially the coupling u of 
the two single QDs to the wire. In this case, the transmission is mainly given by the resonant 
transmission through the wire, and the two single QDs become parts of the reservoirs as it 
is directly seen from ()17|) . 

Moreover, Fig. El (d) shows the transmission through a double QD system consisting of 
two different single QD connected by a wire for the small ratio v/u. Also in this case, the 
transmission probability completely follows the eigenenergies of the double QD. 

All the results obtained for a double QD with one-site single QDs do not show any 
transmission zeros in energy. In order to demonstrate the absence of zeros, we plotted in 
Fig. 121 (a) the energy dependence of the transmission probability for an arbitrary but fixed 
length L of the wire. The model underlying the results of Fig. 01 can be considered as a 
one-dimensional chain of three sites. It is in complete agreement with the consideration by 
Lee that odd and even resonance levels alternate in energy in realistic Id systems so that 
zeros in the transmission probability can not appear. The results of Fig. El correlate also 
with the consideration of a simple two site system jiol . 3|- It has been shown for these 
systems that an architecture of the couplings between system and reservoirs, which violates 
the true one dimensionality of the closed system, gives rise to a zero in the transmission 
probability at a certain energy. 



C. S-matrix for the transmission through double dots: single QDs with many 
states. 

Here, we consider the transmission through a double QD when each single QD of the 
system is presented by two states as shown in Fig. 01 The Hamiltonian (0) of the closed 
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QD 



w 



QD 




□ 




FIG. 4: The double dot system is connected to the reservoirs by the couphng constants v. The 
single dots are coupled to the wire by the coupling constants u. 

double QD consisting of the two single QD and the wire is 

^ ei u ^ 
Q 82 u 
Hb = u u e(L) u u ■ (24) 

U 62 

u El J 

For simplicity we assume that all the coupling constants between the wire and the single 
QD are the same and are given by the constant value u. 

The Hamiltonian ()24|1 is written in the energy representation (jSj), (0}. In order to specify 
the connection between the reservoirs and the single QDs, we have however to know the 
eigenstates of ()24|) also in the site representation. The Hamiltonian of the single QD in the 
site representation is 



V 



/ 



H,= '° . (25) 

The hopping matrix elements Uf, are shown in Fig. lUby thin solid lines. The eigenfunctions 
and eigenvalues are the following 




(jki) = 

We introduce the projection operators 



1 

V2 



I 



1 



P. 



IU(i, 



(26) 



(27) 
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where Bl = 1,2, bn = 1,2, and |1«,) is the one-dimensional eigenstate of the wire. Let Em and 
\m) with m = 1, 5 denote the five eigenenergies and eigenstates of ((211), HB\m) = Em\m). 
The elements of the left coupling matrix are 

{L,E\V\m) = {L,E\VPL\m) = Y,{L,E\V\6b,){6b,\m) 

bL 

= j:iL,E\V\jL){jL\e,,){e,,\m). (28) 

Similar expressions can be derived for the right coupling matrix. Here we used the assump- 
tion that the left reservoir is connected only to the left single QD and the right reservoir 
only to the right single QD. As previously, the reservoirs are assumed to be semi infinite 
one-dimensional wires. Next we have to specify which sites of the left (right) single QD are 
connected to the left (right) reservoir. There are two possibilities. 

(i) Assume the left reservoir is connected only to the first site ji = 1 of the left single 
QD. Then ()28p becomes with account of ()26p 



mE\V\m)=vf-^j:i^,,\m). (29) 

* bL 

A corresponding expression can be written down for the right coupling matrix if the right 
reservoir is connected to the first site of the right single QD. 

(ii) We can assume that the reservoirs are connected to both sites of the single QDs with 
the same coupling constant v. Then the elements of the coupling matrices are the 
following 

{^E\V\m) = vJ^{e^\m) (30) 



provided that the energy level ei is the lowest in energy, see 

The most important difference between the previous d = 1 case and the present d = 2 one 
for the single QD is that the system is now no longer necessarily one dimensional. Therefore, 
zeros in the transmission probability may appear. 

As it was shown above, Eq. (jl(J|) . two eigenvalues of the Hamiltonian Hb coincide with 
the energies ei and £2 of the single QD. The other three eigenvalues of (j^^ can be found by 
solving of a cubic equation. Also the finding of the eigenstates of (j^^ is a formidable task. 
In what follows, we consider therefore the transmission through a system with two states of 
each single QD numerically. 
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FIG. 5: (a) The transmission through a double QD with two identical single QDs that are connected 
by a wire according to Fig. ^ The eigenvalues of Hb are shown by full lines, ei = —1.7, £2 = —1.4 
and e(L) = —1 — L/5 (dashed line), v = 0.3, u = 0.1. (b) The modules of the transmission 
amplitude \t{E,L)\ for the same double QD as in (a) for fixed lengths L = 2.75 (solid line) and 
L = 4 (dashed line). The energies of the two single QDs are shown by circles. The real part (c) 
and imaginary part (d) of the 5 eigenvalues of the effective Hamiltonian as a function of L for 
E = —1.5. Thin solid line: zi, dashed line: Z2, thick solid line: Z3, dotted line: Z4, and dash-dotted 
line: Z5. At L = 2.75 the imaginary part of the third eigenvalue is equal to zero at all energies E. 

In Fig. the transmission probability versus energy E and length L of a double QD 
is shown for the case that both sites of the single QD are connected to the reservoir with 
the coupling matrix elements (jHUj) . The figure shows a zero in the transmission probability, 
indeed, see Fig. El (b). According to Figs. El (c) and (d), the positions and decay widths 
of the eigenstates 2 and 4 of the effective Hamiltonian are independent of the length L of 
the wire while those of the other states depend on L. The state 3, lying in the middle of 
the spectrum, crosses the transmission zero at L = 2.75. Here, the decay width of this state 
approaches zero for all energies E. The transmission zeros will be discussed in detail in the 
next section. Here, we remark only that resonance states with vanishing decay width are 
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considered also by other authors. In [l^, they are called ghost Fano resonances that appear 
in a double quantum dot molecule attached to leads. In atomic physics, the phenomena 
related to resonance states with vanishing decay width are known as population trapping 



221 . They result from the interplay of the direct coupling of the states and their coupling 



via the continuum under the influence of, e.g., a strong laser field. 

IV. ZEROS IN SINGLE-CHANNEL TRANSMISSION 

In Fig. El (a), the transmission through a double QD system with two identical single 
QDs is shown, while Fig. EKb), shows the transmission through one of these single QDs (the 
lower curves correspond to the modules of the transmission amplitude and the upper curves 
to their phases). In the double QD system, the two single QDs are connected to the leads 
and to the wire as shown in Fig. |3 Comparing the two figures, we see that the transmission 
zero of the double QD coincides with that of the single QD. This result is in agreement with 
formula for the single-channel transmission through identical dots. However there is a 
remarkable difference between the zeros in both CcLSGS cLS will be explained in the following. 

Single QD [Fig. ini(b, d)]: The transmission zero of the single QD is due to the destructive 



3,0,0 



and is located between the 



interference of the two neighboring resonance states 
energies of the single QD. It is caused by the unitarity of the 5* matrix with account of 
the fact that the leads are attached to the single QDs which are constituents of the double 



dot system [21|. Around the energy Eq, the transmission amplitude vanishes, t(-E'o) = 0. 
Here Re{t{E)) ~ - Eof while lm{t{E)) ~ ^ - Eq). It holds therefore for 

the modulus of the transmission amplitude \t{E)\ ~ \E ~ Eq\ near Eq, see Fig. El^b), and 
UoT^ 0- Thus, the phase of the transmission amplitude arg{t{E)) jumps by vr at Eq 
according to [12]. The geometrical origin of this phase jump can clearly be seen in Fig. IHfd). 
By pathing through the origin of the coordinates Re(t) = 0, Im(t) = 0, the value arg(t) is 
not defined unambiguously. The phase jumps by ±7r, and the sign of the phase jump is not 
observable. We can call such a zero a first order zero. 

Double QD [Fig. |ni(a, c)]: The zeros in the transmission through the double QD with two 
identical single QDs are of another type. They are of second order since it is ~ |ti(-E)p 

in the vicinity of the energy E = Eq. It follows therefore t{E) [E — EqY near Eq, see Fig. 
El (a). The energy evolution of the real and imaginary parts of the transmission amplitude is 
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FIG. 6: (a) The modules \t{E)\ (lower curve) and the phase arg{t)/TT (upper curve) of the trans- 
mission amplitude for a double QD with two identical single QDs. The parameters are the same 
as in Fig. |S1 and L = 4 as in Fig. [51 (b), dashed line, (b) The modules |ti(i^')| (lower curve) 
and the phase arg{ti)/TT (upper curve) of the transmission amplitude for one of the single QDs 
that is part of the double QD considered in (a). The energy positions of the single QD levels are 
shown by circles at the abscissa, (c) Evolution of imaginary and real parts of the transmission 
amplitude t of the double QD with energy. At the upper right corner, a zoomed fragment of the 
evolution is shown which demonstrates that the evolution has cusp-like behavior in the vicinity of 
the transmission zero, (d) The same as (c) but for the single QD. Here, the evolution is of standard 
type. 

shown in Fig. Efc). In the inset of the figure, the evolution of Re(t), Im(t) at the origin of 
the coordinate system is shown in zoomed resolution. It has a cusp-like behavior, and there 
is no phase jump at all. We present in Fig. IHl (a) the phase behavior of the transmission 
amplitude t that is a combination of two jumps with opposite sign, resulting in a zero phase 
jump at the point E = Eq. This result agrees with the general statement given by Barkay 
et al that phase jumps do not appear when ^^^^^ 1^0= (as in our case) at the energy 
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According to formula (^J, zeros of second order in the single-channel transmission of a 
double QD are given by zeros of first order in the transmission of the single QDs (that 
constitute the double QD) when they are identical. If both single QDs have N energy levels 
then — 1 transmission zeros of second order will appear in the double QD system. A 
numerical computation for the particular case A^ = 5 confirms this conclusion (Fig. [7j). 




FIG. 7: (a) The transmission through a double QD consisting of two identical single QDs with 
d = 5 states that are connected by a wire, versus energy and length of the wire. The eigenvalues 
of Hb are shown by thin Unes while the energy e(L) = 1 — L/2 of the mode in the wire is shown 
by the dashed line, v = 0.5, u = 0.2. (b) Energy dependence of the modules of the transmission 
ampUtude for L = 4. The energy levels of the single QD are shown by circles at the energy axis, 
e L = —2 cos{iTn / Q) , n = l,2,...,5. 

Next we will consider the transmission through a double QD consisting of two different 
single QDs coupled to the wire and to the leads as shown in Fig. |3] When each single QD 
has A^ states, the number of zeros in the transmission through the double QD is, according 
to (0), 2(A^ !)• This conclusion is demonstrated by the results of numerical calculations 
shown in Fig. |H1 Here, the two two-site single QDs are chosen to be different from one 
another what can be achieved by either different coupling constants u between the single 
QDs and the wire or by different energies of the levels of the two single QDs. We have 
chosen the latter possibility. 

In Fig. ISl we see two transmission zeros of first order. The phase jump is — vr at the first 
zero and vr at the second zero. Considering the transition ej, el. with = 1, 2, we see that 
the transmission zeros will approach each other with the consequence that the transmission 
zero turns over into a second order zero. The phase jumps annihilate each other as shown 
in Fig. ini (a). 
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FIG. 8: (a) The transmission through a double QD with different single QDs connected by a wire, 
and (b) the energy dependence of the modules (solid line) and of the phase arg[t)/'K (dashed line) 
of the transmission amplitude for this double QD system at a fixed length L = 4 of the wire. The 
energy levels of the left single QD are = —1.7,-1.4, while those of the right single QD are 
= -1.6,-1.3.. Further, e(L) = -1 - L/5, v = 0.3, u = 0.1. The full lines in (a) are the 
eigenvalues of Hb while the dashed lines are the energies e^, e^, and e{E). 

We consider now the evolution of the modules of the transmission amplitude and the 
corresponding phase shifts when the decay width of one of the states approaches zero. The 
results shown in Fig. El are performed for the double QD system the transmission of which 
is shown in Fig. Altogether with the eigenvalues Zk of the effective Hamiltonian ifefr, Eq. 
fll4|) . as a function of L. The latter ones are related to the poles of the S-matrix. At 
L = 2.75, the third eigenstate crosses the energy of the transmission zero, and its decay 
width Im(z3)/2 approaches zero. As long as L ^ 2.75 and Imlz^) ^ 0, the phase of the 
transmission amplitude varies by vr more or less smoothly, according to the phase shift 
caused by a resonance state with a finite decay width. When L 2.75 and Im(2;3) 0, the 
phase jumps by vr due to the vanishing decay width of the resonance state. Therefore, we 
have also in this case a phase jump of the transmission amplitude by tt. Correspondingly, 
the transmission zero becomes of first order at L = 2.75, see Fig. Efb). That means the 
resonance state whose decay width vanishes when crossing the energy of the transmission 
zero, restores the first order of the transmission zero as well as the phase jump by tt. 

The connection of the two single QDs to the leads and to the wire may differ from that 
shown in Fig. HI When one of the single QDs loses its 2d-character by the manner it is 
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FIG. 9: The energy dependence of the modules (bottom) and of the phase arg{t{E)) /t: 

(top) of the transmission amphtude for L = 2.25, 2.5, 2.75, 3.0, 3.25. The other parameters 
are the same as those in Fig. |S| The transmission zeros are denoted by stars. They are of second 
order. The ordinate is shifted every time by 0.1 when L is changed by 0.25. Ah phases are shifted 
by vr. 

integrated in the double QD system, only transmission zeros of first order appear in the 
double QD. When both single QDs are included as Id dots, the double QD system has no 
transmission zero at all. In this case, the whole system behaves as a Id chain of sites without 
any transmission zeros. 

V. CONCLUDING REMARKS 

From the present study, we conclude that the simple model used by us for the description 
of a double QD system consisting of two single QDs coupled by a wire, is efficient and 
describes qualitatively the features of the resonant transmission through a realistic double 
QD of the same structure. The difference between the transmission through the realistic 
double QD shown in Fig. ^and the model cases shown in Figs. El El El and |H1 consists, above 
all, in the fact that all features of the simple model are multiply repeating in the realistic 
case. These multiple effects in realistic cases are related, obviously, to the large dimension 
of the single QDs and to the higher number of eigenmodes of the wire inside the double QD 
system. 
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The advantage of the simple model described by S matrix theory is that it allows a clear 
discussion of the main features of the transmission, especially of the transmission zeros and 
of the phase jumps related to them. Of course, the model can give only qualitative results 
that do not agree quantitatively with the numerical results (Fig. obtained for the realistic 
double QD system. However there is a large room to further develop the model. First, we 
can take the eigenenergies of the wire for the single channel propagation as 

6„(L) = ^ + ^,n = l,2,... (31) 

where d is the width of the wire. The first term in (j31|) is related to the first-channel 
propagation in the wire. Moreover we can take into account that the eigenstates of the wire 
have the form 

^n(a;) = J|-sin(7rnx/L). (32) 



Then u can be calculated in the same way as v |lfiL \2( 

Un{L) = uoL-^'^i^'^x) U=o,L= UnL-^'\ (33) 

and the Hamiltonian Q reads 

iJ^ = ei(L)|l)(l| + e2(L)|2)(2|. (34) 

In our calculations we have chosen the widths of the wire inside the double QD and the 
widths of leads attached to the QD to be the same. Thereby, the leads and the wire both 
support only single channel transmission through the double QD. Experimentally, however, 
the width of the wire varies, usually, by applying gate voltages |3] while the length of the 
leads remains fixed. In order to model this case, we can consider two energy levels of the 
wire one of which is related to the first channel propagation (j!?T|l and the other one to the 
second channel propagation 

e„.(^) = ^ + -^,n =1,2,.... (35) 

When the wire's width d exceeds the lead's width, we can take one energy level (jHT|) (say, 
n = 2) together with another one from (|35p (say, n' = 1) in such a manner that they cross 
at a certain length Lq, eni^Lo) = e'„(Lo). 

The numerical results obtained for the transmission probability (Fig. ITTHl for such a 
case are the following. For single channel propagation in the wire, the transmission fea- 
tures including the transmission zeros are similar to those discussed above. The number 
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of transmission zeros remains unchanged, but the number of transmission peaks is doubled 
corresponding to the two modes in the wire. These results correspond to those obtained for 
the single channel transmission through a realistic double dot system (Fig. P). The situation 
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FIG. 10: The transmission probability (a,c) and its logarithm (b,d) for a double QD with a wire 
that has two energy levels. The single QDs have eigenenergies = —1.25, —0.75. The coupling 
strengths u between the single QDs and the wire are taken according to l\3'6\i . u = 0.2, v = 0.6. 
In (a) and (b), the first-channel propagation is shown with the two energies of the wire given by 
H31|) . n = 1,2. In (c) and (d), the two channel propagation in the wire is shown with the energies 
ei(L) = 4tt'^/L^ and e2{L) = 1 + 7r^/L^ that are crossing. The dashed lines show the eigenenergies 
of the isolated wire and the single QDs while the solid lines show the eigenenergies of the closed 
double QD system. 

is however another one when the wire can support two channel propagation. In this case, 
the number of transmission zeros is doubled. One of the lines of the transmission zeros does 
not depend on the wire length L (as in the model calculations of Sect. IV), but the other 
one depends on L and crosses the first line of the transmission zeros. The multi channel 
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propagation induced by a variable width of the wire can lead, therefore, to an essentially 
more complicated picture of the transmission zeros. The picture of transmission zeros as a 
whole remains however comprehensible on the basis of the results obtained in the framework 
of the simple model considered above. 

It is experimentally easier to vary the wire's width than its length. It is important 
therefore to mention that the width d and length L of the wire are equivalent for the energy 
levels of the multi-channel wire according to Eqs. (jHT|) and ()35|) . Therefore, the transmission 
probability versus energy and length shown in Fig. can be qualitatively considered as 
the transmission probability versus energy and width of the wire. 
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